!-------------------------------------------------------------------------------
SUBROUTINE reluct(Model, n, X, Y)
!-------------------------------------------------------------------------------
  USE DefUtils
  IMPLICIT NONE
  TYPE(Model_t) :: Model
  INTEGER :: n
  REAL(KIND=dp) :: X(*)
  REAL(KIND=dp), POINTER CONTIG :: Y(:,:)
!-------------------------------------------------------------------------------
  
  Y = 0._dp
  Y(1,1) = 0.001_dp
  Y(2,2) = 1/(sqrt(x(1)**2 + x(3)**2)*1e3 + 1e5)
  Y(3,3) = 0.001_dp
!-------------------------------------------------------------------------------
END SUBROUTINE reluct
!-------------------------------------------------------------------------------
